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In the commented article modifications of tlie Spin- 
Boson Modcl^, for a binary self- regulating gene, have 
been proposed. The new master equations allow two 
different decay rates for free and bound proteins. It 
also presents a vigorous critique of the article 'Self- 
Regulated gene: an exact solution— and 'subsequent 
publications'—"— reducing the finding of exact solutions 
for the spin model to "just a claim" of the authors and 
not as a rigorous result with the strength of a theorem. 

Despite the main criticism on the mentioned solu- 
tions they are really dueling with corrections to the mas- 
ter equations for stochastic gene expression viewed as a 
many-body problem. In the section dedicated to our ar- 
ticle the exactness of our solution is questioned again and 
even a demon, perhaps ironically, is used to explain one of 
the approximations made in the formulation of the Spin 
Boson Model. The task of the demon is the instantaneous 
replacement of the decaying bound protein to the cyto- 
plasm and the inexactness of the solution is a heritage of 
the inexactness of the master equations. The confusion 
between solutions and equations goes until the end of the 
paper when they state that 'our' equations (meaning the 
Spin Boson Formalism) and theirs are dedicated to the 
same problem so one of them must be wrong! Either 
demons or exact equations do not belong to science but 
the concept of an exact solution for a differential equa- 
tion is well defined and exhaustively studied in the field of 
symmetries and integrable systems. Indeed we have been 
able to show that solubility of the Spin Boson models has 
underlying symmetries'"^. 

Both equations contains several approximations and 
are very far from biological details even if we considered 
a simpler induced transcription in E.Coli. The most im- 
portant limitation of this class of models is the bypassing 
of the whole transcription process and even a demon is 
not able to take the information from the DNA, to build 
the mRNA and translate it in the ribosomes instanta- 
neously, putting back the protein to control the promoter 
site. Using the terminology introduced by the authors, 
science has a plenty of demons and the crucial question is 
the selection of proper approximations needed to predict 
and to explain the experiments successfully. 

In section V the authors presents a section named: 



Numerical Validation of The Exact Solution which com- 
pletes the confusion. Exact solutions, by definition, don't 
need a numerical check and are verified analytically. A 
solution can not be validated by a numerical simulation 
or even by an experiment: the occurrence of an agree- 
ment only says that the equations are not wrong. In con- 
trary, the exact analytical solutions are frequently used 
to validate numerical simulations, stability of algorithms 
and so on. 

The authors report the finding of an exact solution for 
the steady-state genetic feedback loop. Examining the 
paper we could not find a complete exact solution for 
the model as announced in the title, but just half of the 
work has been carried out. In fact, they present an exact 
general solution for one probability distribution named 
Gi{z) but they fail to get the other one Go{z) in closed 
form as they recognize explicitly. Mentioning that the 
equation obtained for this component is not of Riemann 
type they only write a first order differential equation 
for it, complaining that "It is difficult to extract an ex- 
plicit solution for Gq{z) by integrating this equation". 
If the authors have read Refsi^i^ carefully, would learn 
that these equations are not Riemann type but from an- 
other family, the Heun equations, and therefore, still in- 
tegrable. 

Consequently the normalization constant is not cal- 
culated obstructing the analyticity of the whole model. 
The underlying reason is that "apparently unknown" in- 
tegrals are needed to obtain the exact solution. They also 
recognize that they cannot calculate the fluctuations as a 
function of the parameters of the model, without the nu- 
merical computation of the normalization, but only the 
fluctuation divided by the mean value, also known as the 
Fano factor. 

The use the generating function technique replaces 
the traditional recursive form of the master equations 
by partial differential equations with the introduction 
of the complex functions Go{z) — J2^=o ^o{n)z^ and 
Gi(^) = S^^Lo where n is the stochastic vari- 

able, the number of protein molecules, and Poin) and 
Pi{n) are the usual probabilities. At steady state limit 
the problem reduces to the solution of the system 

(z - l)ipuGo - G'o) + {9 + auz)Gi - a^zG'o = 0, (1) 
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{z - l){pbGi - G[) - (6 + a,, )Gi+abG'o^O. (2) 
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Gi{z) may be written as a product between Ae''''^^ 
and the KummerM function M(a,/?,w), where 



a = 9 + a,, 



Pu - Pb 



Pu - Pb ~ PbO-b 

+ crf,(cr„ + Pu) 



(1 + Tfc)3 



(3) 

(4) 



{l+(Tb)z-l 
(1 + CTbV 

We obtain Go{z) taking hnear combinations of the Eqs. 
© and ©: 



Go(^) = Ae'"'^^-^) X 



1 + (jf, a 

H 

P« - Pb 



M{a + 1, 13, w) 



M{a,l3,w) 



, (6) 



as we can verify by direct substitution. Normalization 
constant A foUows from probability conservation as usual 
and gives to this model same status of that presented by 
Hornos et.al. 

Throughout their article the authors explain why our 
solutions are not exact and theirs are. The reason is 
that the equations for the spin boson model of Refi^ are 
not exact, some terms are missing and others are "non- 
physical" . Consequently our solution is not exact be- 
cause our equations are not exact. In contra-position 
their stochastic equations are "exact" and therefore the 
partial solution they presented are exact. 

Even though we are not claiming that our solutions 
have the same status, let us consider the exact solution 
for the Schrodinger equation for hydrogen and numerical 
solutions for the Helium atom. Of course the Schrodinger 
equation for an electron in the Coulomb field is not ex- 
act and unphysical if one want to use those terms. The 
Coulomb potential depends only on the radius, relativis- 
tic terms are missing, quantized electromagnetic field is 
absent. The use of the space dependent potential vio- 
lates special relativity and is incompatible with Maxwell 
equations that are Lorentz invariant. An static poten- 
tial between to charges is only possible if we assume the 
existence of a " daemon" which instantaneously tells one 
charge that the other did a small movement. Follow- 
ing their reasoning the hydrogenic solution in terms of 
spherical harmonics and Laguerre functions are not ex- 
act! Of course this is not the case for two reasons:(i) 
there are no exact equations in Science and (ii) we con- 
sider an exact solution for a given equation if we can 
solve it analytically in a closed form. Even if one relax 
the definition of exact solution in closed form to numeri- 
cal evaluation of exact solution, i.e. numerical solutions 
to a target equation, without approximations, calculated 
with desired precision, one cannot call the exact solutions 
presented at Refs.^^^' as just a claim or use quotes to di- 
minish the value of the calculation. In practice the exact 
solutions are powerful because they allow any calculation 



and the study of the parametric equation of one model. 
Furthermore, approximations and physical intuition are 
also easier if we have exact solutions. 

The stochastic equations proposed in their article are 
different from the spin boson ones. Indeed there is no 
parameter that can be used to obtain rigorously one set 
of equations from the other. Both are unphysical and 
approximated. The transcription process, the ribosomes, 
tRNA's and the amplification of the protein number by 
repeated translation of one mRNA have been eliminated 
in both models which is a radical oversimplification of 
any protein synthesis process. Even in the transcription 
of the lac operon of E.Coli controlled by the Lac Re- 
pressor, the comparison with experiments requires two 
stochastic processes'^. In the compared models another 
daemon is required, it carries one mRNA on a hand and 
the ribosomes on the other, translating instantaneously 
and repeatedly the anticodon information. 

The proposed model is welcome and has some inter- 
esting modifications as claimed by the authors. The in- 
stantaneous reposition of protein to the binding mode is 
avoided and the protein degradation when the protein is 
bound is taken into account, even though a new param- 
eter must be included. However, the strict interpreta- 
tion of both models, and, consequently, of the parame- 
ters, is not promising. The transcription and translation 
processes considered here have several time scales. For 
example, the free proteins degradation, the bacterial di- 
vision time, the mRNA degradation time that usually is 
much smaller than the protein degradation, the discrep- 
ancy between mRNA and protein number of molecules by 
more than an order of magnitude in E.Coli^. This means 
that we must have an extra guiding principle to perform 
a physical analysis and the best one are the available ex- 
periments. The critique present here do not diminishes 
the importance of their proposal which can be considered 
an alternative for negative self-regulation. 
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